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Abstract. We consider the classical coupled, combined- field integral equation formulations for 
time-harmonic acoustic scattering by a sound soft bounded obstacle. In recent work, we have proved 
lower and upper bounds on the L^ condition numbers for these formulations, and also on the norms 
of the classical acoustic single- and double-layer potential operators. These bounds to some extent 
make explicit the dependence of condition numbers on the wave number fc, the geometry of the 
, scatterer, and the coupling parameter. For example, with the usual choice of coupling parameter 

^ , they show that, while the condition number grows like k^'^ as A; — >■ oo, when the scatterer is a 

circle or sphere, it can grow as fast as k"^'^ for a class of 'trapping' obstacles. In this paper we 
^\ • prove further bounds, sharpening and extending our previous results. In particular we show that 

there exist trapping obstacles for which the condition numbers grow as fast as exp(7fc), for some 
7 > 0, as A: — > oo through some sequence. This result depends on exponential localisation bounds 
on Laplace eigenfunctions in an ellipse that we prove in the appendix. We also clarify the correct 
^^ , choice of coupling parameter in 2D for low k. In the second part of the paper we focus on the 

) ^ ■ boundary element discretisation of these operators. We discuss the extent to which the bounds 

^_( ' on the continuous operators are also satisfied by their discrete counterparts and, via numerical 

experiments, we provide supporting evidence for some of the theoretical results, both quantitative 
and asymptotic, indicating further which of the upper and lower bounds may be sharper. 



o 






1. Introduction. Consider scattering of a time-harmonic (e~"^* time depen- 
dence) acoustic wave m* by a bounded, sound soft obstacle occupying a compact set 
n C M.'^ {d — 2 01 3) with Lipschitz boundary F, which is such that the complement 
^ . set flf, := M''\ J7 is connected. The medium of propagation, occupying n^, is assumed 

■^ ' to be homogeneous, isotropic and at rest. Under the assumption that u^ is an entire 

C"^ . solution of the Helmholtz (or reduced wave) equation with wavenumber k = uj/c > 

^—^ ' (where c > denotes the speed of sound), we seek the resulting time-harmonic acous- 

tic pressure field u, satisfying the Helmholtz equation 

2. 



en 



^. Am + fc^M == in rie . (1.1) 



O 






This is to be solved subject to the sound soft boundary condition 

u = on r = dne, (1.2) 

and the Sommerfeld radiation condition, which requires that 



S: ^„ifc„«^o(r-('^-i)/2) (1.3) 
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as r :— \x\ -^ cx), uniformly in x := x/r, where u'' :— u — u^ represents the scattered 
part of the field (see e.g. [H]). This problem has exactly one solution under the 
constraint that u and Vm be locally square integrable; see e.g. [26 . 

In this paper we consider the two standard second kind boundary integral equation 
reformulations of (|l.ip - (|1.3p . The first is the indirect formulation 

Ak.vf = 9, (1-4) 

where 

with r; G M\{0} the coupling parameter, I the identity operator and Sk and Dk the 
single- and double-layer potential operators. These are defined for ip G ^^(T) by 



Skifix) := 2 I ^x,y)^iy)ds{y), xeT, (1.5) 



and 



DMx) -.^ 2 J^^^^ ^{y) ds{y), xeT, (1.6) 

with d/dv{y) the derivative in the normal direction, with the unit normal i'{y) directed 
into fie, and $(a;, y) the standard free-space fundamental solution of the Helmholtz 
equation. This is given by 

LH^^\k\x-y\), d = 2, 



47r|a; — y\ ' 



3, 



for x,y € K.'*, x ^ y, where Hq is the Hankel function of the first kind of order zero. 
Finally, g := -2m* jr = 2u''|r. 

The second formulation is the direct formulation 

^'m|^ = /, (1.8) 

where 

A'k,^ := I + D[-iT]Sk, 
with D'/, the integral operator defined, for (p E L'^(r), by 

D'M^) := 2 / ^1^ p{y) ds{y), x G T, 

and 

rill''' 

f{x):=2 — {x)~2ir^u\x), xeT. 
ov 

It is well known (see ^ for details, in particular regarding how classical results can be 
adapted to the general Lipschitz case) that, for ?; ^ 0, Ak^^i and A' k,r) are invertible 
as operators on i^(r), and that 

P'mII-PmII, P'm^'IHPmII- 



(Throughout the paper || • || denotes the L^ norm on F.) 

A question that has received much recent attention in the hterature (see for 
example [21ll[51[Tni[IIl[IIl[IZll2nilII]) is that of determining how the conditioning 
of the two standard integral equation formulations, (|1.4p and (II. 8p . depends on the 
wavenumber k, on the coupling parameter ry, and on the shape of F. Specifically we 
are interested in upper and lower bounds on the (identical) condition numbers of Ak^n 
and A'k,ri, given by 

condA'k,,^ =cond Afc,,, = ||Afc,,,|| ||^fc_J,||, 

and so we are interested in upper and lower bounds on the norms of Ak^n and its 
inverse, and also on the norms ||S'fc|| and ||i?fc||- 

In our recent paper [10], we derived estimates which, to some extent, make ex- 
plicit the dependence of each of these norms on fc, ry and F, with an emphasis on 
understanding conditioning in the important but difficult and relatively neglected 
case where k — > oo. For example, with the usual choice of coupling parameter rj — k, 
while the condition numbers of Ak^rj and A'^. grow like k^'^ as fc — > oo when the scat- 
terer is a circle or sphere [15], we show in [1^ that they grow like fc^/^ for a starlike 
polygon and as fast as fc^/^ for a class of 'trapping' obstacles. In this paper we prove 
further bounds sharpening and clarifying our previous results, in particular studying 
trapping obstacles in much more detail. A main focus of the present paper is also 
the boundary element discretisation of these operators. Our aims here are threefold: 
to provide supporting evidence for some of the theoretical results of [10] and of ^ 
via numerical experiments; to determine how sharp the quantitative upper and lower 
bounds on norms of |10j may be, particularly in the cases where there is a significant 
gap between the two; to determine the extent to which the bounds on the continuous 
operators are also satisfied by their discrete counterparts. 

We begin in Sj2] by summarising the estimates at the continuous level derived 
in [1^, together with previous related results in the literature. Also, in the 2D case, 
we sharpen the estimates from |10j at low frequencies, and prove that the choice of 77 
in [211 ED] (based on analysis for a circular scatterer) guarantees a bounded condition 
number in the limit fc — >■ even for general Lipschitz F. But the main novelty of ^ 
is that we show that there exist trapping obstacles for which the condition numbers 
of Ak^jj and its adjoint grow as fast as exp(7fc), for some 7 > 0, as fc — > 00 through 
some sequence. This result depends on exponential localisation bounds on so-called 
'bouncing-ball' type [19] Laplace eigenfunctions in an ellipse. For completeness we 
provide a self-contained and relatively elementary proof of this exponential localisation 
in the appendix; for eigenfunction localisation results in much more general settings 
proved using related but much more technical arguments see |30| . 

In f|3] we prove results about the relationship between the continuous integral 
operators and their discrete counterparts, i.e. matrices derived from standard Galerkin 
boundary element method (BEM) discretisations. In [j4]we present numerical results 
showing Galerkin BEM approximations to ||Afc_^||, ||^^^||, ||5'fe|| and ||-Dfc|| foravariety 
of obstacles, each for a range of values of k and 77. Finally in ^ we present some 
conclusions. 

The results of the present paper and of [lOj have direct relevance to the numer- 
ical performance of boundary integral methods, since the condition number of the 
discretization of (|1.4p and (|1.8|) appears naturally as a measure of the difficulty of 
computing numerical solutions in practice. Moreover the results in |10| . and more 
particularly our new results on trapping obstacles, have direct relevance to a recent 



detailed fc— explicit numerical analysis of hp boundary integral methods for general 
Helmholtz scattering problems in [53]. There it is shown (for example in [Ml Corol- 
lary 3.18]) that, provided ||^^];|| < Ck^ with C and /3 independent of fc, then an 
hp refinement strategy in which p grows logarithmically in fc and h decreases like 
fc~^logfc yields a Galerkin method which is free from "pollution" (i.e. the error is 
bounded by the best possible error in the finite element space, multiplied by a con- 
stant independent of fc). Our analytical and numerical results are suggestive that 
ll^fc fcll < Ck^ holds for some fc-independent C and /3 not only for starlike obstacles, 
as considered previously in [TT], but also for certain trapping obstacles. But also we 
prove in i j2.5l that there exist Lipschitz obstacles for which the bound ||^^^|| < Ck^ 
does not hold for any C and /?. We also note that the paper [37] (a companion paper 
to [2J) contains new decompositions of the combined potential operators Ak^jj and 
■^'k n "which are crucial in the analysis of the hp methods in |24) . Of key importance 
there is the fact that the decomposition involves certain operators which map into 
spaces of functions which are analytic in a neighbourhood of F. However this analysis 
is rather different in flavour (and has different goals) from that of the present paper. 

We flag that a related and complementary study of the same boundary integral 
equation formulations that we consider in this paper has been carried out recently in 
[5j. That paper includes, similarly to our 21 a numerical study relating to a range of 
geometries of 2D scatterers, but [6] has a different focus, namely an investigation, via 
computation of the numerical range of boundary element discretisations, of conditions 
which ensure that A^^^ is coercive, and how its coercivity constant depends on fc. 

Finally we note that, in a similar vein to our 21 Warnick and Chew [M] |331 132] 
study the conditioning of boundary element discretisations of the single-layer potential 
operator Sk via an approximate theoretical analysis and numerical experiments, ob- 
taining simple explicit approximate upper and lower bounds for the condition number 
as a function of fc and the discretisation step size for several canonical 2D geometries 
(a circle, crack and two parallel cracks) [211 Table 2]. 

2. Bounds on norms and condition numbers at the continuous level. 

2.1. The case of a circle or sphere. Prior to [TU], most research was focussed 
on the case when F is a circle or sphere, in which case Fourier analysis methods are 
possible. 

For the case d = 2, when F is the unit circle, rigorous upper bounds on ||yl/j^^|j 
and ||j4j7„|| for the case t] = k (previously proposed as optimal for conditioning for 
the unit circle when fc > 1 in e.g. [31131 HI]) were derived in [TS] and are that, for all 
sufficiently large fc, 

\\AkM\ < Ck^'\ (2.1) 

II^mII < 1, (2-2) 

with C a constant independent of fc. (Combining (|2.2p with Lemma 12.11 below we 
see that, in fact, ||^^^|j = 1 for all sufficiently large fc.) Although the focus in [T5] 
was on bounding A^^k rather than on bounding the separate components Sk and Dk, 
Lemmas 4.1, 4.9 and 4.10 in [l^ also imply the separate bounds that 

||5fc|| < Cfc-2/3, \\Du\\<Ck^^\ (2.3) 

with C a constant independent of fc. 



For the case d = 3, when F is a sphere of unit radius, it is further shown in |15| 
that, for ah sufficiently large fc, (j2.1[) holds (see also [T7]) and that, for every C > 1, 

for all sufficiently large k. A more refined and fiexible upper bound on A^.r) than (|2.ip 
in the 3D case was recently derived in [5], where it was shown that, for all sufficiently 
large fc. 



\\Dk\\ < C, \\Sk\\ < Ck-^'\ 
for some constant C independent of fc, and hence 

\\Ak,,\\ = \\I + Dk - ivSkW < l + cfl + |r,|fc-2/3 



(2.4) 

(2.5) 

The choice |?7| = k yields the same estimate as (I2.ip . whereas the choice |?7| = fc^/^ 
yields a fc— independent bound for ||A/j_^||. 

2.2. The case of a starlike obstacle. Consider the case when 51 is connected, 
piecewise smooth and starlike, with F Lipschitz and C^ in a neighbourhood of almost 
every a; € F, and 

S- := essinf x ■ vix) > 
xer 

(assuming, without loss of generality, that the origin lies in J7 (0 G J7)). Under these 
assumptions it is shown in [11] that, for 77 € M \ {0}, 



\A 



where 



2 






AS* 



k,ri I 



^ ^ + 1 



<B, 



(2.6) 






2 



d-2 



^*2l 



(1 + 2kRo 
252 r/2 



1/2 



with 



i?o 



sup|a;|, (5+ := esssup a; • z^(a;), S* := essswp\x — {x ■ i'{x))i>{x) 
xer 



xer 



xer 



These assumptions hold, for example, if J7 is a starlike polygon or polyhedron (and 
e 51), and in these cases S^ and S+ are the distances from the origin to the nearest 
and furthest sides of F, respectively. Note that the expression B blows up if fc/JTyj — >■ 00 
or if i5-f/(5_ — > 00, or if 5^\ri\ — > 0, uniformly with respect to the values of other 



variables. If F is a circle or sphere, i.e. F = {x : |a;| = Ro}, then (5_ 



so 



B^Bo 



+ 



1 



^2 



d-2 ^ {l + 2kRaY 



2RW 



1/2 



Rq and 



(2.7) 



In the general case, since 6- < 6+ < Rq and < 6^, < Rq, it holds that B > Bq. 

Based on low frequency asymptotics and numerical calculations for the case when 
F is a circle, it is proposed in [3D] to choose 



n = max ( , k 



(2. 



to minimise the condition number of A^ ,, (and see [3J|3] for some further evidence sup- 
porting this choice). Based on computational experience, Bruno and Kunyansky [71 |B] 
recommend the similar formula that ry = max(6T~^, /c/tt), where T is the diameter of 
ri, on the basis that this choice is found to minimise the number of GMRES iterations 
in an iterative solver. With either of these choices ||^^,J| is bounded uniformly in k 



for fc > for ri starlike. In particular, with the choice (j2.8|) we see that 



\A 



k,ri I 



1 
< - 

- 2 



5_ 



AS* 



[4 



\W + 46i2 



1/2 



2(rf - 2)i?fl 8*^ 



8^0 



1/2 



(2.9) 



where Q = Rq/S-. 

2.3. Upper bounds on | 5fc| , \\Dk\\ and |j^fe.,,| in the general Lipschitz 

case. It is shown in [lOj Theorems 3.3, 3.5, 3.6], under the assumption that the 
scatterer f2 is Lipschitz, that there exist positive constants Ci, i — 1,2,3, dependent 
only on ft, such that 



\\Dk\\ < C^k^'-'y^ + Cs, 
\\AkJ <l + C3 + Csfc^-^-D/^ + Ci|77|fc('^-3)/2, 



(2.10) 
(2.11) 
(2.12) 



for fc > 0. In 2D (d = 2), for the case F simply-connected and smooth, ()2.12p was 
shown previously, for all sufhciently large fc, in [15] . 

Expressions that are in principle computable for the constants Ci, i = 1,2,3, are 
given in |10| . In particular, in the simplest case that F is a straight line of length a, 
the upper bound on |15fc|| is given explicitly by 



ll^fcll <2) 



(2.13) 



The bounds p.l0p - p.l2p are sharp in their dependence on k in the limit fc — ?> 
except for (I2.10p (and so ()2.12p ) in the 2D case. In the 2D case the low frequency 
behaviour is more subtle, as studied previously for the case of a circle in [211 120]. To 
obtain sharper bounds for low k for general Lipschitz F in the 2D case, note that, 
from the power series representations for Yq ^ (9.1.13)], it follows easily that 



Yo{t)~-\ogt 

IT 



for < t < 1, where 7 = 0.577 . 
[H (9.1.18)], this implies that 



2 I, . 1 

<- log2-7 + - 
n 4 



is Euler's constant. Since also | Jo(i)| < 1 for t > 



2i 



4'\t)-^logt 




< 1.03, 0<<<1. 



(2.14) 



r(i) 



— logi -I- 1 -I- 0(1) as t — > O"*", this upper bound is an overestimate 



(Since W^'it) -- 

by not more than 3% for small t.) 

Let So denote the single-layer potential operator in the Laplace case, defined 
by (|1.5p with^{x,y) replaced by <I>o(a;. y) := (l/27r) log(i?o/|2; — yj), for some constant 

6 



i?o > (later we will choose i?o to be some characteristic length scale of F). It is 
a known result (e.g. [33]) that So is a bounded operator on i^(r). Further, (I2.14p 
implies that, where D :— sup^ ,^gp |a; — y\ is the diameter of F, 



$(a;, y) - $o(a;, v) + ^ ^og{kRo) 
for kD < 1. From this inequality it follows that 



<0.26, (2.15) 



11^.11 < ||5o|| + \\S, Soil < \\So\\ + (^n^i^ + 0.26) |r|, 

for kD < 1, where |F| = /p ds is the length of F. Thus, taking i?o = 1 in the above 
result, and combining this bound with ()2.10p . we obtain a refined version of ()2.10|) for 
small k when rf = 2, that 

ll^fell <C7o(l-logfc), forO<fc<l, 

where the positive constant Cq again depends only on fi, which, combined with (j2.1ip . 
gives that 

\\Ak^J<l + C3 + C2k'/'^ + Co\v\{l~logk), 0<k<l. (2.16) 

2.4. Lower bounds on US'/til, ||-Dfc|| and |jAfe^^||. The following lower bounds 
on \\Skl \\Dk\\ and \\Ak,J are derived in [H §4]. 

Lemma 2.1. fT^ Lemma 4-1] In both 2D and 3D, if a part of T is C^ , then 

\\Akj > 1, II^mII ^ 1- 

Theorem 2.2. JKA Theorem 4. 2] In the 2D case, if V contains a straight line 
section of length a, then 

\\Sk\\>\[^+0{k-^) 
V Trfc 



|A,,„||>Hy^-l + 0(HA:-i) 



as fc —> 00, uniformly in r/. 

Theorem 2.3. ]1C^. Theorem 4--4] I"^ the 2D case, ifV is locally C^ in a neigh- 
bourhood of some point x* on the boundary then, for some constants C > and 
fco > 0, it holds for all k > kg and all rj CzM. that 

\\Sk\\ > Ck-^/^ and \\Ak^4 > C\r]\k-^/^. 

More generally, adopt a local coordinate system OX 1X2 with origin at x* and the 
Xi axis in the tangential direction at x* , so that, near x* , F coincides with the curve 
{x* + t*Xi + n*f{Xi) : Xi G R}, for some f G C^{R) with /(O) = /'(O) = 0; here 
t* and n* are the unit tangent and normal vectors at x* . Then if, for some N G N, 
F is locally C^+^ near x* , i.e. f e C^+i(]R), and if also /'(O) = /(2)(0) = • • • = 
f^^'{0) = 0, then there exist C > and fco > such that 

\\Sk\\ > Cfc-(^+i)/(2^+i) and \\AkJ > C\rj\k-<^^+'^/^^^+'^ 

7 



for all k > kQ and all 77 G R. 

In fact, under the conditions of Theorem l2.31 assuming further that /(^+i) (0) 7^ 0, 
we have quantitative lower bounds on \\Sk\\ and ||Afc.^||: 



11^.11 > Cj,{0) fc-(A'+i)/(2^+i) (1 + o(i)), as k 



00, 



and 



/ |r/|Cjv(0)fc-(^+i)/(2^^+i)(l + o(l)), if |7;|fc-(^+i)/(2^+i) ^ 00, 

"'■'■"" - 1 |r/| Cn{0) fc-(^+i)/(2^+i) - ^ + o(l), if |7?| « fc(^+i)/(2A'+i), 



2^2 

as fc -^> cx), where 



2 |/(^+i)(0)| 



1/(27V+1) 



Noting that /"(O) is the curvature at x* , we have the following corollary by 
applying these equations with A^ = 1. 

Corollary 2.4. JIM Corollary 4-5j Suppose (in the 2D case) that T is locally 
C^ in a neighbourhood of some point x* on the boundary and let R be the radius of 
curvature at x* . If R < 00, then, 

\\Sk\\>U^) (2fc)-2/3(i + o(l)), asfc->oo, (2.17) 



and 



lA-.Jl > 



(f)'/' (2^:^/3(1 + 0(1)), ^f\v\k-'/' 



00, 



Mil - I M (I) V3 (2fc)-2/3 _ _^ + „(i)^ ^f 1^1 ^ ^2/3^ 

as fc —> 00. 

We also have the following lower bounds on Iji^fcll- The conditions of Theoreni l2.5l 
are satisfied, for example, if F is a polygon. (Choose x^ to be a corner of the polygon 
and x'^ to be some point on an adjacent side, with F^ a neighbourhood of x^ on the 
adjacent side to x'^ and F^.) 

Theorem 2.5. 110, Theorem 4-6] In the 2D case, suppose x^ and x^ are distinct 
points on T, that F is C^ in one-sided neighbourhoods F^ and F^ of x^ and x"^ , and 
that [x^ — x^) • v(x) = for x € F^ while {x^ — x"^) is not parallel to F^ at x^ . Then, 
for some constants C > and fcp > 0, it holds for all k > k^ that \\D}J\\ > Ck^'^ . 

The conditions of the next theorem are satisfied with A^ = by some pair of 
points x^ and x"^ whenever F is C^. 

Theorem 2.6. '^10, Theorem 4-7] In the 2D case, suppose x^ and x^ are distinct 
points on T, and that, for some N ^ Nq := N U {0}, F is C^ and (7^+^ in one- 
sided neighbourhoods F^ and F^ of x^ and x^ , respectively, and that x^ — x^ is not 
parallel to F^ at x^. Without loss of generality, choose F^ so that, for some e > and 
/eC^+i(R) withf{0)=0, 

F^ = {x^ +tu + f{t)h : 0<t<e} 

where u = (x^ — x'^)/|x"'^ — x^| and h are orthogonal unit vectors, and suppose that, 
for some A^ G No, 

/(")(o) - /(i)(o) = ... = /W(0) = 0. 



Then there exist C > and k^ > such that 

\\Dk\\ > Cfc^/(4^+4) 

for all k > kf). 

2.5. Lower bounds on I^^T J,| for trapping obstacles. In [10] it is shown 
that if 51 is a certain type of trapping obstacle then ||j4j7„|| can be unbounded as 
fc — >■ oo. The type of trapping obstacle considered in ilO^ is an obstacle for which 
there exists points P and Q on the boundary F such that: 

(i) r is C^ in neighbourhoods of P and Q; 

(ii) the line segment joining P and Q lies in f2e and; 

(iii) this line segment is normal to F at P and Q. 
The line segment PQ is an example of a periodic orbit, by which we mean that it is 
the possible locus of a point billiard particle moving in the exterior region ile in a 
straight line at unit speed as on an ideal billiard table, interacting with the boundary 
F according to the usual law of specular reflection (angle of reflection equals angle of 
incidence) . 

The specific class of trapping obstacle discussed in fTO\ is one for which F is a 
straight line locally to both P and Q. Precisely, the following theorem is proved, 
showing that ||^^^|| is unbounded as fc — > oo for some class of trapping obstacles, 
at least provided \ri\ < Ck for some constant C, which is the case of course for the 
standard choice rj — k. 

Theorem 2.7. 110, Theorem 5.1] There exists C > such that, if fie contains a 
square of side length 2a, two parallel sides of which form part of T, and 77 G M \ {0}, 
then 

\\A^!J\ > CkT"(l + l^y\ meN, 

V Km J 

where km '■= rm: /2a. 

Theorem 12.71 relates to the case when the periodic orbit is between straight line 
parts of F. A key idea in its proof is the construction of a quasimode for the Helmholtz 
equation in J7e, by which we mean a function v G H^lTt,,) which satisfies Av + k'^v = g 
with g G L^(f2e) having a small norm relative to that of Av; precisely, the quasimode 
is constructed, dependent on fc, in such a way that if fc = km = rm:/(2a) then 
||ff||L2(a^)/||Av||i,2(f^^') = 0{k~^) as fc — >■ 00. We note that the rate of growth of 
||j4^ II predicted in Theorem 12.71 will be confirmed by numerical calculations in §4, 
and'cf. [U. Fig. 4.7]. 

A periodic orbit between two parallel straight lines is neutrally stable, by which we 
mean that a small initial perturbation in the point billiard's position or direction will 
cause a perturbation to the billiard motion which grows at most linearly with time. 
If the parts of F neighbouring P and Q are curved slightly, so that the periodic orbit 
becomes stable, then the construction of a quasimode becomes possible for which 
ll5llL2(f2^)/||Aw||i2(n^) decreases very rapidly as fc — >■ 00 through some unbounded 
sequence of values (see \T9[ [23] and the references therein) , which leads to a very fast 
growth in H^jT J || as fc — >■ 00 through the same sequence of values. 

We will prove this statement in Theorem 12 .81 below in a case for which a complete 
proof can be given by fairly elementary arguments. This simplest case is that in 
which the parts of F neighbouring P and Q form part of the boundary of an ellipse, 
precisely an ellipse of which PQ is the shortest periodic orbit, in which case the 

9 



quasimode can be constructed by perturbing a so-called bouncing ball mode (see [19| ) 
eigenfunction of the ellipse. This mode can be written down explicitly in terms of 
Mathieu functions and can be shown to be exponentially localised around the stable 
periodic orbit PQ. (For details see the appendix, and for a visualization of several 
of these eigenfunctions see Figure l473l below.) An example of an exterior domain J7e 
and the corresponding scattering object fl which satisfies the conditions of Theorem 
12.81 is the obstacle labeled 'Elliptic cavity' in Figure |4T] below. 

Theorem 2.8. //, for some ai > 02 > 0, fie contains the ellipse E := {(a;i,X2) : 
[xi/aif' + {x2/a2)^ < 1}, and if F coincides with the boundary of this ellipse in 
neighbourhoods of the points (0,±a2), then there exists a sequence < ko < ki < 
k2 < ■ ■ ■, with kjn — >■ 00 as m ^ 00, such that, for some 7 > and C > 0, 

P,-t,J|>Ce^'=- (^1 + M^ , (2.18) 

for T]£E.\ {0} and m^ 0,1, 2,.... 

Proof. In the appendix we focus on a particular subset of the eigenfunctions of the 
Laplace operator with Dirichlet boundary conditions in the ellipse E. These are the 
functions Um.,o G C^iE), m = 0,1,. .., defined by Umfiix) = Mcq (^, (7m)ceo (1^,(7™), 
X € E, where the elliptic coordinates (fi, v) and standard Mathieu function notation 
are as defined in the appendix. The important property of the function Um,o is that it 
satisfies the eigenproblem (jA.ip for wavenumber k = fc,„, where km — 2^Jq,n/{a\ — a\) 
and qm is the (to + l)th positive solution of the equation (|A.6|) in the case n — 0, 
with qm ^ 00 (so that fcm — >■ 00) as m — ?► cxi. It is shown in the appendix (see 
(|A.16P ) that this particular subset of eigenfunctions Wm,o, m = 0, 1, . . ., is a family of 
bouncing ball modes, with Um,Q becoming increasingly localized around the periodic 
orbit O := {(0,X2) : \x2\ < 02} as m — j' 00. 

We will now construct a quasimode Vm on ile by a suitable modification and 
extension oium,o- Let x € C°°(IR^) be compactly supported and such that x{x) = 1 in 
some neighbourhood of O while xi^) = in some neighbourhood of 9£'\F. Abbreviate 
Um,o as Um and define Vm G C^(^e) by Vmix) := x(a;)um(a;), x e E, Vmix) := 0, 
X € il,e\ E. Then, in Vie, 

^^m + kj^V,n = 5m J 

where gm{x) = ior x ^ ile\ E, while 

ffm = Um^X + 2VX ■ VUm 

in E. Let E^ := {x € E : \xi\ > e_} and E+ :— {x € E : \xi\ > e+}, where 
e_ > e+ > arc chosen sufficiently small so that x = 1 in £'\i?_. Then, where ||.||oo 
denotes the usual supremum norm on C{E) and || • II2 the usual norm on L'^{E) and 
II ■ IIl2(-E±) the L^ norm on E±, we see that 

||.9m||2 < ll"m||L2(ij_)||x||co + 2|| Vx|| oo || Vu^ || i2(£;_ ) . 

In the remainder of the proof let C denote a positive constant, whose value does 
not depend on to,, but which is not necessarily the same at each occurrence. By 
(|A.16|) . for some /3 > 0, 

\\um\\mE^)<Ce-f"''-\\um\\2, (2.19) 
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for TO = 0, 1, ... . Further, || Vum||i2(£;_) can be bounded by a constant multiple of 

k,n\\Um\\L^E+), SO that 

hmh < Ck^e-^'^^WUmh, TO -0,1,.... (2.20) 

To see this last claim choose / e C^{M.) such that: (i) 1 > /(s) > for s e M; 
(ii) /(s) — for \s\ < e+; (iii) /(s) = 1 for \s\ > e_; (iv) for some constant M > 0, 
|/'(s)|/(/(s))i/2 < M for aU s eR for which f{s) > 0. (This can be achieved by 
defining / by /(s) = P((|s|-e+)/(e_-e+_)) for e+ < \s\ < e_, where P(t) = tH3-2t).) 
Define X € C^{E) by x(a;) = f{xi), x e E, and note that 1 > x{x) > 0, for a; G E, that 
X(x) = for a; e £■ \ £■+ and xix) = 1 ior x € E_, and that \V x{x)\ / \/ xix) < M 
for all a; e -E for which x{x) > 0. Now, by Green's theorem and since Um is an 
eigenfunction in £', 



lE 



/ Xu^dx ^ - / xu^Au^dx= / V(xu„) • Vu„ da;, 

J_B J_E JB 



SO that 

/ xC^Umfdx < k'^^ ul^dx + M y^ \\/ Um\ Um dx . 

Je J e+ Je 

Applying Cauchy-Schwarz and noting that 2ab < rjo? + r/^^b'^, for all r; > and 
a, 6 > 0, we see that 

for all T] > 0. Choosing 77 — M^^ we see that 

||VM„i||i2(£;_) < Wy^VUmh < C„i||Mm|lL2(_B+), 



where Cm := y/M'^ + 2fc^ 
Next note that 



\Vm\\l> y-lidx =\\Um\\j-\\u^\\l2(E^)> -\\Urn\\l, (2.21) 



for all sufficiently large to,, by (|2.19p . Thus, and noting ()2.20p . 

II Au„j2 = ll^mWrn - ffmlU > fcTO||Wm||2 - llffmlb > ^||Wm||2, 

for all sufficiently large m. Combining this bound with (I2.20p we see that 

ll5m||L2(o.) _ ||,g™||2 <cfc;;ie-^'='"<Ce-^'='", TO = 0,1,..., (2.22) 



for some < 7 < /3. 

To see that (|2.22l) induces exponential growth of ||A^^ ||, we proceed as in the 
proof of dni Theorem 5.1] and define «;, e C{R'^) D H^^^iM.'^) by 

vl{x) := / $fe„ (a;, y)g„,iy) dy, x e R^ (2.23) 

Je 
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where $fe^ denotes the fundamental solution $ of the Helmholtz equation in 2D 
in the case k = k^. Then we can view !;,„ e C^{ilf^) as the total field for the 
problem of scattering by the obstacle Vt in the case when v]^ is the incident field. For 
defining u^ := Vm — W™ it holds that Aw^ + k'^vf^ = in fig, that vf^ satisfies the 
Sommerfeld radiation condition (since u^ does and Vm is compactly supported), and 
that w^ = — u^ on. F. It follows, arguing as in the proof of jTOl Theorem 5.1], that 



a; 



dVr, 



^m jV ' 



du 



Jrr 



(cf. (HH)), where 



f„,(x):=2^{x)-2irjv:^{x), xeT. 



(2.24) 



Since H^^,^^^ „|| = ll(^fe,„ ,,) ^11 7 our proof of (|2.18p will be completed if we can show 
that, for some constant 7 > 0, 



dvn 



dv 



L^T) 



>Ce^'^-[l + ^] ||/™||i2(r), 



(2.25) 



for m = 0, 1, . . . and rj ^M.\ {0}. 

To see that (|2.22p implies (|2.25p . we use (|2.23p . and we also apply Green's repre- 
sentation theorem [TS] to Vm to give that 

v„,{x)^ I ^k^{x,y)grn{y)dy+ I ^k„A^,y)^iy)ds{y), x e il^. (2.26) 

Using the bound (e.g. fTUl) that \H'^^\t)\ < ^2/(7ri), for t > 0, which implies that 

|*fc„(a;,y)|<(87rfc„|x-y|)-i/2, 
we easily deduce from (|2.26p that 



\\v,n\\2<Ck-y'i\\g^h + 



dv„ 



di 



LHF) 



for TO = 0, 1, ... . Combining this bound with (|2.20l) and (|2.2ip we see that 

dv„ 



dv 



L^r) 



> CkK^u. 



m\\2, 



(2.27) 



(2.28) 



for TO = 0, 1, ... . Similarly, it follows from p.23p that 

\\vln\\mr)<Ck;„'/^9rnh, 

and that 

Vvl{x)^wl^\x)+wi^\x), 
where, for x G M^, 

u^';n\x):^ / S/^<^o{x,y)g.m{y)dy, w't^^x) := / \7^ {^k^{x,y) - <^o{x,y)) g„^{y) dy 

J E J E 
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(l/27r) log(l/|x — y\) is the standard fundamental solution 
Now, from standard mapping properties of Newtonian po- 



and(cf. §2.3)$o(a;,2/):= 
of the Laplace equation. 

tentials, it holds that Wm^ € H^{E), with ||wm''||ffi(_E) < C||(7m||2- Hence, by the 
boundedness of the standard trace operator from H^{E) to H^^^{dE) D L^{dE), it 
v^{dE) < C'||(7m||2- Further it holds (see e.g. p^ equation (3.9)]) 



follows that 
that 



\w, 



(0), 



|V,($o(x,2/)-<i>fe,„(a;,y))|<C. 



\x~V\ 



from which it follows (cf. (|2.28p ) that 

Hence 

dvl 



dh 



L\T) 



< CkU^9n.h, 



and combining this bound with (I2.28|) and the definition (j2.24p of /,„ , we see that 






Il5m||2- 



Finally, combining this bound with (|2.20p and p.27p . we see that 



I/™ II 



L2(r) 



< Ckme 



'Pk„ 



dv„ 



dv 



LHT) 



which implies that (I2.25[) holds for 7 < /3. D 



2.6. Choice of ry for low k. Although the main focus of this paper is on condi- 
tioning in the limit as fc — > 00, for completeness we briefly address the limit fc — >■ in 
this section. Conditioning in this limit was explored carefully already in the papers 
[2TI I20] where, for the case when F is a sphere or circle, precise asymptotic calcula- 
tions were made of the choice of rj which minimises cond Af^jj in the limit A: — >■ 0. The 
recommendations in these papers are for a circle/sphere of unit radius, and imply for 
a circle/sphere of radius i?o > that the optimal choices of 77 are 



77 



^^+Oik^logk), d-3, 

{7r2 + 4(log(A:/2) + 7)'}"'/'{l + 0{k^ log fci?o)}, d = 2, 



(2.29) 



where 7 = 0.577 ... is Euler's constant. We will explain in this section why these 
choices, for any Rq > 0, ensure a bounded condition number of Aj^-jj as fc — >■ in the 
case of general Lipschitz F. 

To understand this limit we need to recall what is known about integral equation 
formulations for the Laplace case fc = 0. Let $0 denote the fundamental solution of 
the Laplace equation, given simply by ()1.7|) with fc = in the 3D case, and defined as 
in §2.31 in the 2D case. Let Sq and Dq denote the single and double-layer potentials 
in the Laplace case, defined by equations (jl.Sp and (|1.6p with S, D, and <& replaced 
by 5*0, Dq and $o- It is a fairly straightforward calculation (see e.g. [lU] for the detail 
in the case of Lipschitz F) that 



\Dk~Do\\^ 



Oand \\Sk-So\ 
13 







(2.30) 



as fc — ?► in the 3D case, and that the first of these resuhs holds also in the 2D case. 
In the 2D case the limiting behaviour of Sk is more subtle. We see from (|2.15p that 



Sk~So + ^ log(fci?o) T 



^0 (2.31) 

as A; — > where T is the finite-rank integral operator defined by 

Tcj)ix) = f (l>{v)ds{y), X e r. 



The following limiting behaviour of Ak^^i is clear from (|2.30p and (|2.3ip . 
Lemma 2.9. As k ^ Q, 

Afe,„ = / + Do-ir;S'o(l + o(l)) 

in 3D, while 

Ak;„ ^I + Do + iv^ log(fci?o) T - b]So{l + o(l)) 

ZTT 

in 2D. Thus, unless 

"^^X 0((logfc)-i), d = 2, "^/^^U, 

it holds that \\Ak^n\\ — > oo as fc — )■ 0. On the other hand, if, for some Cq G M, 

7? ^ Co ask^O, (2.32) 



in the case d — Z or 



in the case d = 2, then 



T]— log(fci?o) -^ Co as k^O, (2.33) 

ZTT 



||^fc,j7 - ^oll ^ as k ^ 0, 
where 

. ^ ( I + Dq - icoS'o, d = 3, 
° ■ ^ I + D0 + icoT, d = 2. 

The above lemma, coupled with the following theorem, makes clear that it is 
appropriate to choose r/ for low k so as to satisfy (j2.32p or (12.33^ . for d = 3, 2, choosing 
Co 7^ in each case. This choice of 77 ensures that cond A/j,, remains bounded in the 
limit as A: — >■ 0. Clearly, one such choice of 77 is (|2.29l) . 

Theorem 2.10. Where Aq is as defined in Lemma \2.9[ it holds that, for co 7^ 0, 
Aq is invertible as an operator on H'^{T) for < s < 1, in particular as an operator 
on L^(r), while Aq is not invertible for co = 0. Thus, if (J2.32I) or (j2.33p hold in the 
cases d = 3 and d = 2, respectively, then, as fc — > 0, H^jT^H = 0{1), if cq ^ 0, while 

Proof. The last sentence follows immediately from standard operator perturbation 
results and (|2.33p and (I2.32p once the first sentence is proved. In the case co = it 
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is well known that Aq is not injective, having a non-trivial null space which is the 
set of constant functions, see e.g. [HI Theorem 6.20], [3Tj. To show invertibility of 
Aq for Co 7^ we note first that, by interpolation, it is enough to show invertibility 
on H''{T) for s = and 1 [2^. Further, since the difference Aq — Ak.ri is a compact 
operator on i^(r) and on H^{T) (see e.g. the proof of Theorem 2.7 in [9]) and since 
Ak^n is invertible, it holds that Aq is Fredholm of index zero on i^(r) and on H^{T)^ 
so that it is invertible if and only if it is injective. Moreover, since Aq is Fredholm 
with the same index on H^iT) and L'^iV), and L^iT) is dense in 7?^(r), it follows 
from a standard result on Fredholm operators (see e.g. [551 §1]); that the null-space 
of Aq is a subset of H^iT) D H^/'^iV). In the case that F is C^ that there are no 
non-trivial functions in the null-space of Aq in C{T) is shown in [35J Theorem 6.24] 
in the case d = 2 and in [T^ Theorem 3.33] in the case d = 3. In the case when F 
is Lipschitz the same arguments can be used to prove injectivity of Aq in if ^/^(F), 
replacing the mapping properties of layer potentials in classical function spaces in 
[121 [22] with those in Sobolev spaces in [26] (cf. the proof of Theorem 2.5 in [9 ). D 



2.7. Bounds on condition numbers and choice of ?/. In this section we 
bring together the results from the sections above and explore their implications for 
the conditioning of Ak^n, and what this then implies regarding the choice of rj to 
minimise cond A^ ,,. We have already noted in ii2.2l and !j2.6l recommendations made 
in the literature regarding the choice of rj, mainly based on study of the case when F 
is a circle or sphere. Overwhelmingly (see e.g. [IHliUliSl 111 (IZllZliSl US ) the guidance 
is to take ry proportional to k for all but small values of k. The choice of r/ for small 
k has been discussed already in ^M One choice of ry, recommended by Kress [20] 
for the 3D case, that we have studied in i j2.21 is 77 = max(l/(2i?o), fc). This choice, 
by Lemma [2.91 above, is not suitable in the 2D case for low k, since with this choice 
IJAfc^^lJ — >■ cx) as fc — )■ 0. An alternative choice, which satisfies (|2.33p with cq y^ 0, and 
which we will use for computations in 21 is 

^_ r (Eo(l-log(fci?o))-\ 0<fci?o<l, (2.34) 

[^ k, fci?o ^ 1- 

Here Rq is a length scale of the scatterer ft; we choose Rq as defined in i j2.2l in 21 
The following theorem, which follows from (HU), ^^, ([232]) . (jTTe]) . and Theorem 
12.101 is a sharpening of results in [TOl Section 6] . 

Theorem 2.11. Suppose that F is piecewise C^ and starlike, in the sense of ^2.^ 
and that rj is given by p.8|) in the case d — 5, by ()2.34p in the case d = 2. Then, for 
some constant C > 1, 

1<\\A^',\\<C, l<||A,,,||<C(l + fc('^-i)/2), 



that 



condAfc,^<C2(l + fc('^-i)/2) 



for all k > 0. In the case d ~ 2 we have a sharper lower bound for k large, so that, 
for some c > 1, 



c-\l + fci/3) < cond Ak,rj < c(l + fci/2), (2.35) 



for k > 0. 
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For the case of a circle or sphere we saw in W2.1\ that the above upper bounds are 
not sharp; with the proposed choices of rj the sharper bound holds that 

cond Ak^rj < C(l + fc^/^). 

We will investigate, in the 2D case, which of the bounds in (j2.35p is sharp in 21 We 
will also investigate the alternative choice for r/ for large k proposed in [5] , namely to 
take \r]\ = fc^/^. It follows from p.5|) and (j2.6|) that, when F is a sphere, this choice 
of rj also implies 

cond Afe,^ < Cfc^/^ forfc>l. 

We will explore whether this estimate holds for 2D geometries in 21 Note that for 
starlike polygons and the choices of rj indicated in Theorem 12.111 it follows from 
Theorem 12.21 that, for some c > 1, 

c~\l + fc^/2) < cond Ak^^ < c(l + fc^/2), 

for fc > 0, i.e. it is the upper bound in (|2.35p that is sharp in this case. We will 
illustrate this in the numerical results in 21 

For trapping obstacles, in the sense defined in ^12.51 faster rates of growth of 
cond Ak^ri are inevitable. The following result is deduced in |10| . by combining The- 
orems [^2] and [2Jl 

Theorem 2.12. \10[ equation (6.13)] Suppose that the conditions of Theorem 
are satisfied. Then, for some C > 0, where km is as defined in Theorem\2.7\ 



cond^,„.,>Cfcr" (l + H(l + fc„0-V2) (l + M^ 



m e N. 



Thus, if rj — c(l + k^^), for some constants c and p, then, for some constant C > 0, 
cond Afe,„,^ > (7(1 + fc^), meN, 

with q = 9/10 forO<p< 1/2, q = p + 4/10, for 1/2 < p < 1, and q = 14/10 for 
p>l. 

For trapping obstacles satisfying the conditions of Theorem l2.8l the situation with 
regard to conditioning is much worse: the condition number must grow exponentially 
as k increases through some sequence of wavenumbers. 

Theorem 2.13. Suppose that the conditions of Theorem \2. 8\ are satisfied. Then 
there exists a sequence < fco < fci < ^2 < . . ., with k„i — > oo as to — > cxd, such that, 
for some 7 > and C > 0, 

cond Afc„,,, > Ce'^"^, 

forr]eR\ {0} a^rf to = 0, 1, 2, . . . . 

Proof. By Lemma [^TT] and Theorems 12. 31 and l2.8[ there exists a sequence < fco < 
fci < fc2 < . . ., with km — > 00 as TO — >■ 00, such that, for some 7 > and C > 0, 

cond Afc„,„ > Cmin (l, Mk-^^) e^'^" [l + ;^ 

But for \rj\ < k„i this implies that cond^fc^.^ > ^Ce^^"^ while for jryj > km this 
implies that cond Ak^^n > ^Ckm e''''^™, and the result follows. D 
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3. Discrete level. In this section we explore the relationship between ||Afc ,J| 
and 11^^^, II and the norms of discrete versions of these operators, specifically the 
norms of matrices arising from Galerkin discretisations. 

Let Xm C i^(r) be a finite-dimensional subspace with Pjv : ^'^(r) -> X^ the 
corresponding orthogonal projection. Let V^ be a bounded linear operator on L^{V)- 
Then, given y £ L^(r), a Galerkin method for solving the equation 

Vx ~ y 

for X £ i^(r), is to seek xn € Xn such that 

PnVxn = PnV- (3.1) 

Let {01, . . . , (/)jv} be an orthonornial basis of Xj^, define V/y '■ Xj^ ^^ Xj^ by Vn '■= 
PNV\xr,, and let T^ : X^ ^ C^ be defined by 

TnX = [(x, 0i) • • • (x, <i)N)T ■ 

Then Tjy is an isomorphism, indeed an isometric isomorphism if we give C^ the 
standard Euclidean norm || • ||2. Further (I3.ip is equivalent to 

V^TnXn = TnPnV, 



where 



V" := TnVnT^^ 



is the linear operator on C whose matrix representation (that we denote also by 
V^) is the Galerkin matrix V^ = [(^0^, 00]- Clearly 

\\Vn\\ = WW (3.2) 

(where we use || • || on the right hand side to denote the matrix norm induced by the 
vector norm || • II2), since both T^ and T^^ are isometrics. Also Vat is invertible if 
and only if V^ is invertible and, if they are both invertible, then 

ll^^'ll=ll(0-'ll- 

Now we need to determine the relationship between ||VAr|| and ||V^||. We first require 
the following result. 

Lemma 3.1. If W is a hounded linear operator on L'^iV) and Pi,P2,... is a 
sequence of orthogonal projection operators with Pn4> ~^ fof o,ll 4> € -^^(r); then 

\\PnWPn\\^\\W\\, 

as N ^ 00. 

Proof. Let A = lim inf Ar„j.oo IIPjvW^-PatII and choose a monotonic increasing se- 
quence Ni,N2,. . . of natural numbers with WPn^^Pn^ || — > A as A; — > cx). Then, for 
every (f) £ L^(r), 

\\W(b\\ = lim \\PN,WPN,(b\\ < lim ||Pjv.W^-PjvJ|||<^|| - A||0||, 

and hence ||W^|| < A = lim inf jv-s-oo 1 1 ^w W^^w 1 1 • On the other hand, we have 

limsup||PjvT^fW|| < ||VF||, 
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since \\PnWPn\\ < ||^7v||||W^||||^w|| = ||W^|| for every iV, and hence 

\\W\\ < liminf llPArTyPArll < limsup||PArTyPAr|| < \\W\\. 

Thus hniAr_j.oo llP/vTyP/vll exists and is equal to \\W\\. D 

Clearly HP/vV^PtvH = l|K/v||i so that it follows from (13.21) and Lemma [3. II that 

WVW^WVr.W^WPNVP^W^WVW, (3.3) 

as TV — >■ oo. 

In the case that V = I + C with C compact, it holds moreover that ||(V^^)^^|| ^ 
\\V~-^\\ as iV — > oo if y is invertible. To see this, note first that, by Lemma [3.11 if V 
is invertible, 

\\PnV-'Pm\\^\\V-'1 
a,s N —i' oo. Next, let Vn = I + PnC and note that 

Vn — Vn\xni 

and that, since P/v converges strongly (i.e. pointwise) to the identity and C is compact, 
PnC converges in norm to C, so that H^— VatH — >■ 0. It follows from standard operator 
convergence results that Vn is invertible for all sufficiently large TV. But then it follows 
that also Vat = VatIxn is invertible (as an operator on Xn)- Indeed, injectivity of 
Vn is clear by injectivity of Vn- To see surjectivity, take -0 € Xn and note that, 
by surjectivity of Vn, there exists a e i^(r) with -0 = VN<i> = + PnCcJ), so 
that (p = ip — PnCcJ) e Xn and hence, VniJ) — V-"- This argument also shows that 
y-i = V^Vjv- Further, \\V^'^ - V'^^W ^ as iV ^ oo, so that 

\\PnV-^Pn ~ PnV^^PnW = \\Pn{V-^ - V^^)Pn\\ 

<\\Pn\\\\V-^~V^'\\\\Pn\\ 

= \\v-'-v^'\\^Q, 

as N —i' CX3. Hence, as A^ — >■ oo, 

ll(^^)-'ll = WVn'W = WPnV^'PnW = \\PNV-'PN\\+oil) ^ \\V-^. (3.4) 



Equation (j3.3[) applies to the Galerkin boundary element method discretisation 
of all the operators we have discussed in the previous sections, in particular to Sk, 
Dk, and A^ ,,, provided that the sequence of approximation spaces Xn is chosen so 
that P/v converges pointwise to the identity. It is enough to check that this pointwise 
convergence holds on some dense subset, for example to check that 

\\Pn(I) - 011 = inf II0JV - 011 ^ as iV ^ oo, for every e C{T). 

Equation p.4p applies to the operator Afe,,, if F is C"'^, for then A^,, has the form 
Ak,r) = I + Ck.-q with Cfc,,, compact 16J. For general Lipschitz F, it is not known 
whether p.4p holds, indeed it is not even known for any Galerkin method that V^ is 
invertible for all sufficiently large N. 
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4. Numerical results. In this section we compute ||V^^|| for V — Sk,Di^,Ai. 



and II (V^^) ^11 for V = A^^i, each for a variety of obstacles, and we compare the 



computed values with the upper and lower bounds on the corresponding continuous 
operators as described in ^ The aim is to provide supporting evidence for some of 
the theoretical results described in [J2j both quantitative and asymptotic, and to give 
some indication of which of the upper and lower bounds may be sharper, particularly 
when there is a significant gap between them. We also seek an indication of the 
extent to which the bounds on the continuous operators are satisfied by their discrete 
counterparts. 

We present results for rj = k for all geometries under consideration, and we also 
present results for 77 = k'^^^ for certain specific examples. As we have discussed in 
W2.71 the choice ry = fc is widespread in the literature, e.g. [31 IH [TSl [T71 [5^, and this 
choice is supported by our own preceding analysis. The interesting choice ij = k^'^, 
proposed in [S], is also supported by some of the above analysis; for example we have 
seen in ij2.7l that, for a spherical scatterer, cond Ak^jj increases at the same rate as 
fc — > 00 whether 77 is proportional to k or proportional to fc^'^. 

Although our main focus is on larger values of fc, for two examples we also inves- 
tigate the limit fc — > 0, presenting results for rj = k and for 77 given by (|2.34p . 

In each example the boundary F is piecewise C°° , that is F = U^i r^-'-' with F'^-'^ 
a C°° arc. We denote the length of F^^^ by Lj, and divide each F'-'^ into Nj segments 
T^', i = 1,. . . ,Nj of equal length iFj-^'^l = Lj/Nj. We then define the orthonormal 
basis functions by 

^' ^ ' 1 0, otherwise, ^^---^^h- 

Now whilst in theory \\V\\ — limjv->.oo ll^^lli in practice we can only compute 
||T^^|| for a finite value of N. In order to justify the assumption that our choice of 
Nj is sufficiently large we fix Nj oc k, choosing the constant of proportionality on the 
basis of some simple model experiments. In particular, for the case that il is a circle 
the eigenvalues of A^ ,, are known explicitly, with corresponding formulae in terms 
of the eigenvalues for ||Afc_^|| and ||^^|,|| (see [lOl §2] for details). Thus for a circle 
we can compare our computed approximations to ||^/c,r;|| and ||^^J,|| with the known 
values. For this example we found that ten basis functions per wavelength gives a 
relative error of approximately 1%, and thus in each example we choose Nj « ^^ ^ . 

We present numerical results for the obstacles shown in Figure HTT] More detailed 
descriptions of the obstacles are provided below. For each obstacle and for each 
operator 14- we also compute the algebraic growth rate p under the assumption that 
||Vfc|| = CkP, for some constant C > 0. Assuming this formula holds, we can estimate 
the value of p from two successive values || T4- || and || Va:^_|_i || by 

^l llVfcll j 

p= v;^\ ■ (4-1) 

log (^) 

In the cases where p « 0, p is not shown in the tables. In all of the estimates detailed 
below, C and Cj, j = 1,2, . . ., denote unspecified constants independent of k and 77. 

4.1. Circle. For our first example, we consider the unit circle. From (|2.17p 
and (|2.3I) we know that, for fc > 1 

(32^)-l/3fc-2/3(l + 0(1)) < ll^fcll < Cfc-2/3^ 
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Circle 



Ellipse 



Kite 





Crack 



Square 



Thin Rectangle 



Rectangular cavity 



Elliptic Cavity 




Fig. 4.1. Obstacles corresponding to numerical experiments. 



where here and throughout this section o(l) denotes a term which vanishes in the 
hmit as k ^ oo. Theorem [^^ (lO)) and (12.111) imply that, for fc > 0, 

whilst we know that the sharper upper bound (|2.4I) holds in the case of a sphere, 
that IJDfclJ < C. The numerical results in Table H?T1 for the corresponding boundary 
element matrices suggest that this sharper result, proved for a sphere, appears to 
be applicable for a circle as well; we observe for the discrete approximations that 
\\Sk\\ ~ k~^/^ and \\Dk\\ ~ fc° {^ in this section indicates that the ratio of the left 
hand side to the right hand side is approximately constant in the limit fc — > oo). The 
quantitative lower bound on H^fcH from (|2.17l) is clearly a lower bound in Table WA\ 
underestimating the true norm by a factor of about 6.5. 

From Lemma [^m Corollary [Ol (P^ and (P?T^ we know that, for fc > 1, 



f- 



1/3 



) (l + o(l)) <\\Ak^k\\ < C^k^'\ 
l<Pfe,fc2/3||<C2A:i/2. 



Note that Corollary 
Il^fc,fc2/3||, as fc ^ oo 



gives the lower bound (327r)-i/3 - 1/(2V2) w -0.139 on 
which is clearly less sharp than the lower bound on \\Ap. J.2/3II 
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k 


{32TT)-y-'k-''r-i 


\\Sk\\ 


P 


\\Dk\\ 


5 


7.355x10^^ 


5.240x10-1 




1.144 


10 


4.633x10^2 


3.152x10-1 


-0.73 


1.114 


20 


2.919x10^2 


1.997x10-1 


-0.66 


1.084 


40 


1.839x10-2 


1.246x10-1 


-0.68 


1.079 


80 


1.158x10-2 


7.798x10-2 


-0.68 


1.076 


160 


7.297x10-3 


4.884x10-2 


-0.68 


1.075 


320 


4.597x10-3 


3.076x10-2 


-0.67 


1.072 


640 


2.896x10-3 


1.935x10-2 


-0.67 


1.071 



Table 4.1 
Circle. Norms of Galerkin BEM approximations to Sk and Di^, and p values given by 1)4. Ij l. 



from Lemma [^H The upper bound on || A^^2/3 1| for the case of a sphere is, from (12. 5p . 
||^fc,fc2/3|| < C3. The numerical results in Table |4?2] suggest that this sharper result 
also holds for a circle; the results suggest ||Afc^fc2/3|| ~ fc°, and that ||^fc,fc|| ~ k^^^ 
as expected. The quantitative lower bound on ||Afe.A;|| from Corollary 12.41 is a lower 
bound in Table 1121 underestimating the true norm by a factor of about 7. 



k 


V327r/ 


Uk.kW 


P 


u-dw 


Il^fe,fc2/3|| 


IM-i II 


P 


Bq,,2/3 


5 


0.37 


2.663 




0.986 


2.016 


0.995 




3.82 


10 


0.46 


3.233 


0.28 


0.987 


1.993 


1.056 


0.09 


4.49 


20 


0.58 


4.021 


0.32 


0.987 


1.981 


1.260 


0.26 


5.38 


40 


0.74 


5.030 


0.32 


0.987 


2.000 


1.701 


0.43 


6.56 


80 


0.93 


6.271 


0.32 


0.987 


1.999 


2.039 


0.26 


8.06 


160 


1.17 


7.859 


0.33 


0.987 


1.990 


2.694 


0.40 


9.98 


320 


1.47 


9.883 


0.33 


0.987 


1.998 


3.407 


0.34 


12.40 


640 


1.85 


12.419 


0.33 


0.987 


2.000 


4.307 


0.34 


15.49 



Table 4.2 
Circle. Galerkin BEM approximations to \\Af^^^\\ 



and \\A 



fc,77 1 



By Lemma I^H ||A^.^|| > 1, which combined with ^^ implies that \\A 



k,k\ 



1 



for all k sufficiently large, and the numerical results in Table |4~21 show this behaviour. 
The bound for general starlike obstacles applied to the circle, i.e. (|2.7|) . gives that 



I^mII ^ I 



^ I fc^ I (1 + 2fc) 

jy2 27^2 
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1/2 



=:B, 



0,r;- 



Note that B, 



0,fe 



2.5 (in fact it holds that 2.5 < i?o,fc < 2.6 for the range of k in Table 
14.21) . and that Bq ^,2/3 ~ -\/3 fc^^^ as fc — t' 00. We see from Table IT2l that i^o,?; appears 
to be an upper bound for the discretisation of ||^^„|| as predicted, overestimatimating 
by a factor of about 2.5 for the larger values of k when 77 = fc, by a factor of about 
3.6 when r] = k"^'^ . 

We note from Table |4?2] that, for this example, the condition number cond Ak^n = 
IJAfc^^ll ll^^^ll appears to be slightly numerically smaller for 77 = fc2/3 than for -q = k. 
It appears that, for both choices of 77, cond A^^n increases approximately in proportion 
to fci'3j though this is less clear in the case ry = k"^'^. 

4.2. Ellipse. Next we consider the ellipse given by (2cosi, isint), t G [0,27r]. 
The more specific results of t ^2.1l do not apply in this case, and for upper bounds 
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on ||5fc|| and ||-Dfc|| we have only the results for general Lipschitz T of H2.3I The 
inequalities (|2.17p and (12.101) imply that, for fc > 1, 

(4^)-i/3fc-2/3(i + o(l)) < \\Sk\\ < Ck-^'"^, 

the lower bound larger than for the case of the circle as the maximum radius of 
curvature {R — 8) is larger. Theorem 12.61 and (12. lip with A^ = imply that, for 
fc>0, 

Ci<\\Dk\\<C2k^''' + C^. 



Inspecting the numerical results in Table 14731 we see that the quantitative lower bound 
on ||S'fe|| from (|2.17p is clearly a lower bound for the norm of the discretised operator, 
underestimating the true norm by a factor of about 6 at the highest wavenumbers (cf. 
the results for the circle). The numerical results for ||Dfe|| suggest that ||-Dfc|| ^ kP ^ i-C- 
that the lower bound on ||-Dj.|| is sharp, while it appears from the numerical results 
that llS'fell -fcP, forpw -0.6. 



k 


(4^)-l/3fc-2/3 ||5fc|| p 


ll^^ll P 


5 

10 

20 

40 

80 

160 

320 

640 


1.471x10-1 6.692x10-1 
9.267x10-2 4.143x10-1 -0.69 
5.838x10-2 2.730x10-1 -0.60 
3.678x10-2 1.803x10-1 -0.60 
2.317x10-2 1.209x10-1 -0.58 
1.459x10-2 8.029x10-2 -0.59 
9.194x10-3 5.269x10-2 -0.61 
5.792x10-3 3.427x10-2 -0.62 


1.458 

1.591 0.13 
1.671 0.07 
1.760 0.08 
1.819 0.05 
1.877 0.05 
1.919 0.03 
1.942 0.02 



Table 4.3 
Ellipse. Galerkin BEM approximations to \\Sk\\ 



and \\Dh 



Now turning to Table 14. 4[ note that Lemma 12.11 Corollary 
that, for fc > 1, 



and (prr^ imply 



1/3 

f-] (1 + 0(1)) <||A,,,|| 



< 



1< IIA 



fe,fe2/3| 



Clfcl/2, 
< C2k^/\ 



(4.2) 



(Note that Corollary [121 gives the lower bound (47r)-i/3 - 1/(2^2) « 0.077 on 
II Aj. j.2/3 II , as fc — >■ cx). So in (14.21) we have used the sharper estimate from Lemma [2Tl ) 
The numerical results in Table 14.41 suggest that ||Afe,fe|| ^ k'P for p sa 0.4, that 
||j4j.j.2/3|| ~ fc°, and that the lower bound in (|4.2p is an underestimate by a factor 
approximately 2.5. 

Lemma I^TT] and (|2.6p imply that, for fc > 1, 



1< IIA 



fe.fel 



<C7i, 1<\\A 



-1 

fc,fc2/3 



< C2k^'\ 



(4.3) 



The numerical results in Table W^ suggest that ||Aj, ^|| w 1 for fc > 5. The values of p 
corresponding to ||yl 



fc,fc2/3l 



are rather variable, but the average of the last six values 
for p is 0.34, approximately consistent with the upper bound C2fci'3. 

As for the case of the circle the condition number cond A^^^ — \\Ak^n\\ ll^fcill 
appears to be numerically smaller for 77 — fc2'3 than for rj — fc, by nearly a factor 2 
for the higher values of fc in Table 14.41 
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k 


ii^r 


ll^fc,fc|| 


P 


ll^.:lll 


Uk^kv4 


ll/l-i II 

Il"^fe,fc2/3ll 


P 


5 


0.736 


3.507 




0.987 


2.417 


0.996 




10 


0.927 


4.267 


0.28 


0.987 


2.473 


1.024 


0.04 


20 


1.168 


5.589 


0.39 


0.987 


2.554 


1.300 


0.34 


40 


1.471 


7.317 


0.39 


0.987 


2.599 


1.662 


0.35 


80 


1.853 


9.751 


0.41 


0.987 


2.580 


1.986 


0.26 


160 


2.335 


12.902 


0.40 


0.987 


2.582 


2.548 


0.36 


320 


2.942 


16.906 


0.39 


0.987 


2.634 


3.387 


0.41 


640 


3.707 


21.972 


0.38 


0.987 


2.689 


4.287 


0.34 



Table 4.4 
Ellipse. Galerkin BEM approximations to ||Aj; ,j|| 



and \\A 



fc,77 I 



1.5 


r— ^ 


/ 
/ 
/ 


1 


^ \ 


A\ 


0.5 


\ 







/ 
/ 
/ 
/ 




-0.5 


/ 


y 


-1 


r 


^^ 


-1.5 


/ , 





-0.5 0.5 



Fig. 4.2. The 'kite' shape and notation of ^4-3\ 



and ||Afc^2/3l 



4.3. Kite. Estimates for \\Sk\\, \\Dk\\, \\Ak,k\\, W^k.kW, \\Ak,kV^\ 
for the 'kite' shape given by (cost+0.65cos2i — 0.65, 1.5 sini), t g [0, 27r], are shown in 
Tables l4?5] and l476l for various k. Our upper bounds from Sj2]that apply in this case are 
identical in their dependence on k to those for the elhpse. However, in contrast to the 
circle and the ellipse, our lower bounds imply significant growth in ||Dfe|| and || A;. ^.2/3 1| 
as k increases. Specifically, applying Theorem l2.3l to the non-convex kite, with N — 2 
and X* the point of infiection on F labelled x^ in Figure IT2l and recalling (|2.10p and 
(Hm, it follows that, for fc > 1, 



Cik-y' < \\Su\\ < C2k-^'\ 
Cik^'^ < \\AkM\ < C2k^'\ 

Clfcl/15<||A,fe2/3|l<C2fcV2. 



Moreover, applying Theorem 12.61 to the non-convex kite, with N = 2, x'^ one of the 
points of inflection on F, and x^ the other point on F intersected by the tangent at 
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x^ (see Figure 



and recalling (j2.1ip . we have, for fc > 1, 



The numerical results in Tables H3] and W^ provide some support for these estimates. 
The lower bound on H^fcH seems sharper than the upper bound, although in fact the 
behaviour of \\Sk\\ appears to be rather similar to that for the ellipse. On the other 



\A 



k,k'^/^\ 



hand, the behaviour of ||l?fe|| and 

ellipse, with, approximately, \\Dk\\ '^ k^^^ and ||Aj._;.2/3|| ^ fc^/^ 
Since the kite shape is starlike, satisfying the assumptions of 



is very different from that seen for the 



AM hold on II A 



k,k\ 



and II A 



fc,fc2/3l 



the same bounds 
as for the case of the ellipse. We observe simi- 
lar behaviour to that of the ellipse, namely that H^^rj^H ~ 1 and that the bound 
Il^fcfc2/3ll < C2fc^/^ appears sharp. In Tableli^ cond A,, ,,2/3 ~ fc^, with p « 0.47, 
a faster rate of growth than cond Ak^k ~ fc^/^, so that, for the larger values of k, 
cond Ak,k < cond Aj, ^,2/3. 



k 


11 ^^11 


P 


\\Dk\\ 


P 


5 


6.591x10-^ 




1.810 




10 


4.365x10-1 


-0.59 


2.169 


0.27 


20 


2.758x10-1 


-0.63 


2.686 


0.29 


40 


1.712x10-1 


-0.65 


3.160 


0.27 


80 


1.074x10-1 


-0.65 


3.616 


0.25 


160 


6.759x10-2 


-0.66 


4.160 


0.24 


320 


4.400x10-2 


-0.65 


4.760 


0.24 


640 


2.866x10-2 


-0.65 


5.437 


0.23 



Table 4.5 
Kite. Galerkin BEM approximations to \\Sk\\ and \\Dk\\ 



k.k\ 



^ 1 



\A 



k,k\ 



P 



\A 



\A 



fc,fc2/3| 



P 



\A 



P 



5 

10 

20 

40 

80 

160 

320 

640 



3.720 

4.766 

6.151 

7.513 

9.316 

11.563 

14.337 

18.387 



0.36 
0.37 
0.29 
0.31 
0.31 
0.31 
0.36 



0.987 
0.987 
0.987 
0.987 
0.987 
0.987 
0.987 
0.987 



2.764 
3.243 
3.602 
4.041 
4.394 
4.891 
5.414 
6.030 



0.23 
0.15 
0.17 
0.12 
0.16 
0.15 
0.16 



1.029 
1.027 
1.177 
1.511 
1.883 
2.354 
2.954 
3.671 




0.20 
0.36 
0.32 
0.32 
0.33 
0.31 



Table 4.6 
Kite. Galerkin BEM approximations to \\Ai^^k\\, \\A^ f_\\, \\Af_ 



fc2/3 



and \\A 



-1 I 

fe,fc2/3l 



4.4. Crack. This numerical example is distinct from the others in that F is an 
open arc, the straight line from (0,0) to (0,1), and the only theory which applies 
from ^ are the upper and lower bounds for \\Sk\\- Computations of the Galerkin 
BEM approximations to ||5fc|| are shown in Table 14771 for various k. From (|2.13p and 
Theorem O it follows that 



^ + 0(fc-^)<||5.||<-|=. 
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(4.4) 



The results in Table H771 clearly demonstrate that ||S'a;|| ^ fc^^/^ (of. [32]), a slower rate 
of decay than for the circle, ellipse or kite, and that the values of \\Sk\\ are bracketed 
between the quantitative upper and lower bounds in (J4.4I) , nearly coinciding with the 
lower bound values. 



k 


ll^fell 


P 


v/l/(7rfc) 


2^1/{7rk) 


5 


2.649x10-1 




2.523x10-1 


5.046x10-1 


10 


1.817x10-1 


-0.54 


1.784x10-1 


3.568x10-1 


20 


1.266x10-1 


-0.52 


1.262x10-1 


2.523x10-1 


40 


8.960x10-2 


-0.50 


8.921x10-2 


1.784x10-1 


80 


6.326x10-2 


-0.52 


6.308x10-2 


1.262x10-1 


160 


4.472x10-2 


-0.50 


4.460x10-2 


8.921x10-2 


320 


3.162x10-2 


-0.50 


3.154x10-2 


6.308x10-2 


640 


2.236x10-2 


-0.50 


2.230x10-2 


4.460x10-2 



Table 4.7 
Crack. Galerkin BEM approximations for \\Sf^\\, and theoretical lower and upper hounds. 



4.5. Square. Computed estimates for ||5fc||, H-DfcH, ||^, 



kM 



\A-A\ 



\A 



k,k^/'>\ 



and \\A 



-1 I 



for the square of side length two are shown in Tables 14.81 and 14.91 



below, for various k. The theoretical upper bounds from g2]that apply are identical 
to those for the ellipse and kite examples above. In particular, since the square is 
starlike satisfying the conditions of §2.2) the bounds (|4.3|) apply. However, as F now 
contains a straight line segment we have different lower bounds on ||S'fc||, ||^fc,r/|l and 
\\Dk\\ compared to the ellipse and the kite. Specifically, applying Theorem 12.21 and 
recalling (|2T0| and ((2T2l) . it follows that, for fc > 1, 



-A;i/6 



\l^ + o{k-') <\\s,\\< ck-''\ 

1-1-0(1) <\\Ak,k\\< Ck^'\ 
l + 0(fc-i/3)<p^.^^.,,3|, <Cfci/2. 



(4.5) 
(4.6) 
(4.7) 

(4.8) 

It appears from Table HT51 that \\Sk\\ ~ fc-i/2^ as expected, and the quantitative lower 
bound in (|4.5p appears to be sharp, underestimating |jS'fc|| by only 3% at the highest 
frequency. It also seems that ||I?fe|| ^ k^/^, indicating that the lower bound in (|4.8p 
is sharp in its dependence on k. 

The results in Table |4!9] suggest that ||j4i;.fc|| ~ fc^, with p w 1/2, as expected 
from (|4.6p . It appears that \\A-^ k'^/A\ i^ increasing roughly like k^/^ . The quantitative 
lower bound in (14.71) is seen to be a lower bound for the Galerkin BEM discretisation 



Applying Theorem 12. 51 and recalling (|2.1ip . we also have 

Cifci/4 < \\Dk\\ < C2fci/2. 



of ||Aj,j,2/3|| in Table 221 underestimating ||A^,j.2/3|| by about a factor 3.5 at the 
highest frequency. As in the cases of the circle, ellipse, and kite, II^^T],!! ~ 1 for all fc, 
while 11^4-, 2/3 II is increasing as fc increases, though the rate of increase is somewhat 
erratic (from (|2.6p . we recall that ||^~!_2/3ll ^ Ck^^^). However, it appears that 
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k 


y7k 


ll^fcll 


P 


\\Dk\\ 


P 


5 


3.568x10-^ 


5.784x10-1 




1.316 




10 


2.523x10-1 


3.353x10-1 


-0.79 


1.488 


0.18 


20 


1.784x10-1 


2.137x10-1 


-0.65 


1.730 


0.22 


40 


1.262x10-1 


1.428x10-1 


-0.58 


2.018 


0.22 


80 


8.921x10-2 


9.760x10-2 


-0.55 


2.389 


0.24 


160 


6.308x10-2 


6.723x10-2 


-0.54 


2.825 


0.24 


320 


4.460x10-2 


4.667x10-2 


-0.53 


3.346 


0.25 


640 


3.154x10-2 


3.259x10-2 


-0.52 


3.972 


0.25 



Table 4.8 
Square. Galerkin BEM approximations for \\Sk\\ 



and \\Di. 



cond Af. ^.2/3 ~ fcP, with p ~ 0.6, a faster rate of growth than cond A^.k ~ k^^"^ , and, 
for the largest value of fc, cond A^^k ~ cond Af.f.2/3. 

Where _Ro > is some length scale of the scatterer, it follows from Theorem 12. Ill 
that choosing 



r,^f1* -.^l/iRoil-HkRo)), 



ensures ||v4fc^,J| and ||^j^„|| remain bounded as A: — >■ 0. This is true for any Lipschitz T 
and (rather arbitrarily) we choose this example to illustrate this numerically. Define 
Rq as in W2.'2l so Rq — -\/2 for this particular scatterer (taking the origin at the centre 
of the square). With this choice of Rq, we show in Table HTTUl norm computations 
We see that, while ||^ri. || seems to blow up for small values 



for small values of k. 

of fc, the values of 11^^^1*11 remain essentially constant, and cond A^,,,* appears to 
be approaching a limit of about 3.7 as k — >■ 0. For the computations in Table l4T0l 
deviating from the element sizes used in the other calculations, we discretised each 
boundary line of the square with 500 equal length elements so as to be able to resolve 
the blow up in the norm of the inverse of Ak.k with some accuracy. 



k 


\\AkA 


P 


ll/l-i|| 


v/l-^ 


1/6 _ I 


^4fc^fc2/3 


P 


IM-i II 


P 


5 


3.089 




1.024 




0.043 


2.237 




1.116 




10 


3.611 


0.23 


1.024 




0.171 


2.375 


0.09 


1.172 


0.07 


20 


4.608 


0.35 


1.023 




0.315 


2.536 


0.10 


1.280 


0.13 


40 


6.032 


0.39 


1.023 




0.476 


2.865 


0.18 


1.484 


0.21 


80 


8.117 


0.43 


1.023 




0.656 


3.187 


0.15 


1.904 


0.36 


160 


11.068 


0.45 


1.023 




0.859 


3.600 


0.18 


2.315 


0.28 


320 


15.253 


0.46 


1.023 




1.087 


4.131 


0.20 


3.065 


0.41 


640 


21.177 


0.47 


1.023 




1.342 


4.777 


0.21 


4.064 


0.41 



Table 4.9 
Square. Galerkin BEM approximations for \\Ak^k\\> ll^fcfcll) ll^fcfc^/sl 



and \\A 



-1 I 

fc,fc2/3l 



4.6. Thin rectangle. The bound (|2.6I) suggests that, even if the obstacle J7 is 
starlike, ll^^i |j may blow up as the obstacle becomes very thin when 5+/5- is small. 
In this example we investigate the extent to which such a blow up happens. Where 
5 is the thickness of the scatterer, it is reasonable to expect this to be a problem 
when k5 <^ 1, for if K{x,y), x,y ^ F, denotes the kernel of the integral operator 
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k 


Pfe,/c|| 


Pmi 


P 


Uk.v'W 


u^:,4 


10"^ 


1.608 


3684.85 




1.738 


2.106 


10-4 


1.608 


465.66 


-0.89 


1.723 


2.105 


10-3 


1.608 


61.63 


-0.88 


1.703 


2.104 


10-^ 


1.608 


9.04 


-0.83 


1.680 


2.100 


10-1 


1.612 


2.11 


-0.63 


1.693 


2.081 


1 


2.391 


1.88 


-0.05 


2.534 


1.871 



Table 4.10 
Square. Galerkin BEM approximations for ||ylfc fe||, \\A 



1 II 
fc.fcll 



ll^fc,^* 



and \\A 



k,rj* 



Dk — vrjSk and x± £T are adjacent points on opposite sides of a thin part of F, then 
K{x+,y) « K{x-,y), y € r, so that the integral operator should be badly conditioned. 
To explore the extent to which this is a problem, and the extent to which the bound 
(|2.6I) reflects actual behaviour in this limit, we show estimates for ||Afe_fe|| and H^^ ^H 
for a rectangle with side lengths 2 and 0.02 in Table 14.111 below, for various small 
values of k. 

Using the notation of i)2.21 for this example we have S* = 5^ = 2, (5_ = 0.02, and 
Ro = V4.0004, and ^^ tells us that 



I^mI 



<B, 



where B is as defined in ^ 32. 21 However, we see from Table ICTI that this bound is a 
gross overestimate, at least provided we choose rj carefully. For the definition of B 
implies that, whatever the choice of 77 G M\ {0}, 



1 
^>2 



(5_ 



45* 



(5_ 



g*2 



> 2-^ = 2 X 10^ 



for this geometry. If we choose rj 

blowup as k 

similar). But, if we choose rj — rj 

[231] that P^J,| 

growth in Table gUl but H^-^.l 



0, as for the case of the square (indeed the values of ||A^^ 



k for small k then, indeed, we see significant 

are 
l/Ro{l — ln(fci?o))j we know from Theorem 
must stay bounded as fc — )■ 0. In fact we see some mild, logarithmic 
is never larger than 26 for the range of k shown. 
For the computations in Table 14.111 we used 100 elements on each boundary 
segment. Comparing the algebraic rates p associated with ||^^jj|| for the first three 
wavenumbers k — 10"^, 10"^, 10"'^ it appears possible that the discretisation does not 
fully resolve ||^^^|| for k = IQ-^ and that the exact value may be higher. However, 
due to convergence issues of the underlying singular value decomposition for the norm 
computation no finer discretisation could be used here. 

4.7. Rectangular cavity. In the last two numerical examples we explore trap- 
ping domains, as studied theoretically in i j2.5l The rectangular cavity in Figure 14.11 
is defined by the polygon with the following coordinates: po — 0, pi = (— c, 0), 

P2 = i-c,-£), P3 = {£,-£), Pi = {£,2c-£), P5 = {-c,2c-£), pe = (-c,2a), 
Pj — (0,2a). Here, a — tt/IO, c = 1, and £ — c ~ a. The width of the cavity is 
2a = 7r/5. Hence we expect resonance values in the negative half of the complex 
plane close to fc = 5m, n G N (see [61 Figure 5.12] for numerical computations of exte- 
rior resonances for this cavity). Since for k — 5m the width of the cavity is an integer 
multiple of half a wavelength. Theorem 12 . 71 applies and implies that, for k = 5m and 
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k 


Pfc.fell 


IIV^II 


P 


II^M-II 


11^.-.;. II 


P 


10"^ 


1.978 


3087.76 




1.978 


25.165 




10-4 


1.978 


924.25 


-0.53 


1.978 


21.533 


-0.07 


10-3 


1.978 


128.19 


-0.86 


1.978 


17.581 


-0.09 


10-^ 


1.978 


45.28 


-0.45 


1.978 


14.078 


-0.10 


10-1 


1.978 


14.43 


-0.50 


1.978 


10.384 


-0.13 


1 


2.050 


4.53 


-0.50 


3.310 


3.570 


-0.46 



Table 4.11 
Rectangle of side lengths 2 and 0.02. Galerkin BEM approximations to WA^; fe|| 

\\Af, „. II , and \\A'^^ 



-1 

fc,TJ* 



m e N, 



i^^,;ii > ck'"'. 



oo, 



\A 



k,k\ 



^fc,fc2/3l 



will grow, a behaviour not 
for various k are shown in 



for both rj ^ k and r/ = fc^/^. Hence, as k 
observed in the previous examples. 

Estimates for \\Ak,k\\, Pfc,lll, Il^fc,fc2/3|| and \\A 
Table \4A2\ below. As expected, the value of ||^^];|| grows at a rate at least as fast as 
fc^/i°, and cond Ak^k grows at approximately the rate fc^/^ predicted by Theorem l2.12l 
The growth of ||^r 1.2/3 II 9'iid cond Aj^ 1^2/3 are more erratic. Theorem 12.51 implies that 



fc2/3| 



> 



\\Dk\\ > Ck'^/^, and |15'fe|| < Ck^/^ by (l2lTO . so that, for this obstacle, P^ 
Cfci/4. Thus the theoretical lower bound for cond Af. f.2/3 is cond Af. f.2/3 > Cfc^'^/^°, 
a slower rate of growth than that proved for cond Ak^k- Nevertheless, cond Ak,k < 
cond Af. f.2/3 for the larger values of k in Table l4T2l 



k 


UkM\ 


P 


II4-MI 


P 


WAk^kV^W 


P 


Pfc.Wall 


P 


5 


4.835 




1.969 




3.575 




3.075 




10 


5.201 


0.11 


3.121 


0.66 


3.580 


0.00 


5.990 


0.96 


20 


5.629 


0.11 


5.539 


0.83 


3.594 


0.01 


9.348 


0.64 


40 


6.182 


0.14 


10.322 


0.90 


3.790 


0.08 


19.029 


1.03 


80 


8.112 


0.39 


19.774 


0.94 


4.223 


0.16 


28.528 


0.58 


160 


11.066 


0.45 


38.351 


0.96 


4.788 


0.18 


173.563 


2.60 


320 


15.254 


0.46 


75.156 


0.97 


5.483 


0.20 


277.480 


0.68 



Rectangular cavity. 

i-^ II 



Table 4.12 
Galerkin BEM approximations to 



^fe,ft2/3| 



and 



4.8. Elliptic cavity. The elliptic cavity in Figure |4~T1 is defined by two elliptic 
arcs. The first one is parameterised as (cost, 0.5 sin t), t e [—(j)o,(j)o], with 0o = 
77r/10 and the second arc is defined by (1.3cosi, 0.6sin<), t e [— (/)i,^i], with 0i = 
arccos(Yi3 coscjjo). As discussed in ^2.5[ we expect large values of H^^ ^H at k values 
corresponding to so-called 'bouncing ball' eigenmodes of the inner ellipse, which has 
semi-axes 1 and 0.5. Some of these modes, discussed in detail in the appendix, are 
shown in Figure 14.31 As k grows, these modes become more localized around the 
centre of the ellipse, i.e. around the stable periodic orbit (see >j2.5l and the appendix). 
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k =9.9771201566136298 




k, „=28.807002784875433 

4,0 




k„ =60.218097688523919 




k,, =91.632551202864647 

14,0 




Fig. 4.3. Some exponentially localised modes of the ellipse associated with zeros of the radial 
Mathieu function of order 0. In the notation of the appendix, the modes plotted are Um,o for 
m = 1,4,9, and 14, corresponding to the wavenumbers km,o, "* = 1,4,9, 14. 



Estimates for ||Afejt|| and H^jT],!! for the four wavenumbers from Figure l4?3l are 
shown in Table 14.131 We expect, from Theorems 12.81 and 12.131 exponential growth 
of ll^^^ll and cond Ak.k with k. Comparing the values of 11^/7^11 for the two low- 
est wavenumbers a large growth is clearly visible. However, the growth of ||^fcj.|| 
then levels off. It may be that the given discretisation with about 10 elements per 
wavelength is not sufficient to resolve the large norm of the inverse, or that, due to 
discretisation error, the resonance wavenumber k is shifted somewhat from the theo- 
retically predicted value. Repeating the computation with approximately 20 elements 
per wavelength, we see the computed values for H^^ j.|| for larger k are significantly 
larger now. However, the growth of the inverse still levels off, and we do not ob- 
serve exponential growth even with this finer discretisation. We note, however, that 
with this finer discretisation we see, at the highest value of fc, a condition number 
cond Ak,k ~ 140, 000, hugely larger than the values observed in any of the previous 
examples, including the rectangular cavity. 



k 


UkM\ 


P 


ll/l-MI 


Aj^^ J, (refined) 


9.977 


5.595 




67.3 


70.1 


28.807 


7.294 


0.25 


2829.9 


9182.7 


60.218 


8.820 


0.26 


5265.0 


13373.8 


91.633 


10.144 


0.33 


6543.8 


14258.6 



Table 4.13 
Elliptic cavity. Galerkin BEM approximations to ||j4(; fe|| and ||A^j.||. For \\A'^^\\ the results 
of a refined computation with twice the number of degrees of freedom and approximately 20 elements 
per wavelength are shown. 
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5. Conclusions. In this paper we liave, in fj^l summarised what is known re- 
garding upper and lower bounds on the norms of the acoustic single- and double-layer 
potential operators, Sk and Dk, and the combined layer potential operator Ak^n^ with 
an emphasis on how these bounds behave as a function of frequency, and the influ- 
ence of the shape of the boundary. We have also proved sharper upper bounds on 
ll^fcll and IJAfc^^ll for low k, have summarised what upper and lower bounds on ||^^„|| 
are known, and have shown that exponential growth of |1^^„|| is possible as fc — > oo 
through some sequence of wave numbers, in the case of a certain class of 2D trapping 
obstacles. Finally, we have discussed the condition number cond ^fc,??, proving that it 
remains bounded as fc — > with appropriate choices of the coupling parameter 77, and 
showing that, while it increases as fc ^- 00 only as fast as fc^ " for a circle or sphere, 
and at the rate fc^" for a starlike polygon, it grows exponentially, as k increases 
through some sequence, for certain trapping obstacles. 

In fJ3] we have explored the implications of these results for Galerkin BEM dis- 
cretisations of these operators, showing that the norms of the Galerkin BEM matrices 
converge to the norms of the operators that they discretise, as the mesh is refined, 
and provided an orthonormal basis is used. Convergence to ||^jr„|| of the norm of the 
inverse of the matrix corresponding to A^jj has also been proved in the case that T 
is C^. Thus we expect that the norm bounds at the continuous level in §2 will apply 
also at the discrete level if the mesh is sufficiently refined. 

This has been confirmed in 21 where we have explored a range of numerical 
examples, including shapes that are convex (both smooth and non-smooth), non- 
convex but starlike, and non-starlike trapping obstacles. The quantitative upper and 
lower bounds stated in fj2]are found to be upper and lower bounds also at the discrete 
level, and to be rather sharp in many of the examples. In almost all cases the observed 
rate of growth of norms and condition numbers as k increases is in accordance with 
the possible range of behaviour suggested by the upper and lower bounds from fj2l 
with the rates of growth mainly closer to the lower bounds of ^12. 41 The exception 
is that the exponential growth predicted in Theorems 12.81 and 12.131 is not observed 
numerically in the 'Elliptic cavity' example, at least at the discretisations we use. On 
the other hand, the condition numbers observed in this case are as high as 140,000 
and the values of ||^j7Jj|| as high as 14,000, which contrasts with 11^4^7^11 ~ 1 in all the 
cases where the scatterer il is starlike. 
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Appendix: Eigenmodes of the Ellipse. In this appendix we summarise k;ey 
properties of eigenmodes of the Laplacian in an elliptical domain. Suppose ai > 02 > 
and let E — {(xi,X2) € R^ : (xi/ai)^ + {xija-if' < 1} be the ellipse with semi-major 
axis fli and semi- minor axis 02. Then we study in this appendix certain Laplace 
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eigenmodes u G C^{E) satisfying 

Aw + fc^u = in £;, u = on dE, (A.l) 

for some fc > 0. 

Let a = \/ a{ — d\ — aie, where e = \J\ — a^/a\ is the eccentricity of the ehipse, 
and introduce elhptical coordinates (/i, J^), defined by 

Xi = a cosh /i cos i^ and X2 — a sinh /i sin z^, 

in terms of which 

E — {{a cosh /i cos v, a sinh /i sin i^):0</z</^0: 0<i^<27r}, 

where jjlq := tanh^^(a2/ai). It is well known (see e.g. [35]) that the Laplace operator 
separates in elliptical coordinates, and that in this coordinate system the Helmholtz 
equation can be written as 

Jf^ + ^fT + ^^^^(sinh^ M + sin^ i^)]u = 0. (A.2) 

Seeking separation of variables solutions in the form u{x) = M(y^)iV(i/), we see that 
(jA.2p implies that N satisfies the circumferential (or standard) Mathieu equation 

A^"(i') + (a-2gcos2i^)iV(i^) = 0, (A.3) 

while M satisfies the radial (or modified) Mathieu equation 

M"(/x) - (a - 2qcoah2fi)M{fi) = 0. (A.4) 



In these equations 



^{kar = \{ka,ni-{a2/a^r) 



and a is a separation constant. The solutions to (jA.3[) that we are interested in are the 
solutions of period 27r, satisfying iV(0) = Af(27r) and 7V'(0) = 7V'(27r). Since N{-v) 
is such a solution if N{iy) is, it is clear that we may restrict attention to periodic 
solutions of (JA.3P that are either even or odd. 

For this paper it is enough to focus on the even periodic solutions, so that we seek 
solutions N of (|A.3p which satisfy A^'(O) — N'{tt). This is an eigenvalue problem in 
which the eigenvalue is the separation constant a. Standard Sturm-Liouville theory 
tells us that, for each value of the parameter g > 0, there are a countable number 
of eigenvalues a — a„((7), n — 0, 1, ... , with ao(g) < ai{q) < . . . and a„(g) — >■ oo as 
n — >■ cx) (we use in this appendix the standard notation for these eigenvalues and the 
corresponding eigenfunctions, see e.g. 14') §28.2(v)]). Further, to the eigenvalue a„(q) 
there corresponds a unique (to within multiplication by a constant), real- valued eigen- 
function N{i'), which is usually denoted by ce„(i/, q) (for the standard normalization 
of ce„ see ^4j §28.2(vi)]). The standard Sturm-Liouville theory tells us that cen(i^, q) 
has precisely n zeros in (0,7r). It is easy to see that ce„(^, g) ± ce„(7r — i/,q) is also 
an eigenfunction corresponding to a = ««(?); thus the uniqueness of the normalised 
eigenfunction implies that, for i/ € K and n > 0, 

ce„ (j^, <z) — ce„ (tt — i/, q) if n is even, ce„ (;/, q) = — ce„ [tt — i>,q) if n is odd. (A. 5) 
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Given that a — an{q), for some ri > 0, and that N — ce„ (•,(/), it is a standard 
resuh (e.g. [33]) that u{x) = M{^)N{v) satisfies the Helmholtz equation (jA.l[) in E if 
and only if M G C^ (R) is an even function that satisfies (IA.4I) . This uniquely specifies 
M to within multiplication by a constant. The standard notation for this (real- valued) 
solution is Af(/i) — Mc;j^'(/i, g); see jTH §28.20(iv)] for the standard normalization. 
Thus we see that u{x) — M{^)N{v) = McJ^'^' (/i, g)cc„(zy, q) satisfies the full eigenvalue 
problem (lA.ll) if and only if 

McW(mo,<z) = 0. (A.6) 

The complication in computing eigenmodes of the ellipse (for methods see |35[[29] ) 
is that it is a multi-parameter spectral problem: to satisfy (|A.6p we have to find a 
pair (a, q) such that, simultaneously, (|A.3|) has a periodic solution and (|A.4p has a 
solution which is even if N is even and which vanishes at jjlq. Neves [2^ gives a proof 
based on multi-parameter spectral theory that for each pair (tti, n) € {0, 1, . . .}^ there 
exists a unique qm,n > such that (jA.6|) holds with Mcjj {■,qm,n) having in zeros in 
(0,/io)- The function 

u{x) = Um,nix) := Mcll\n, qm,n)cen{v, qm,n) (A.7) 

is then an eigenfunction of (jA.l|) for k = km,n '■— \/Mm,n/0'- It is well known (e.g. [12] ) 
that the eigenvalues of the Laplace operator have infinity as the only accumulation 
point, so that fcm,n — > oo as m + n -^ oo. 
For some i/q G (0, 7r/2) let 

Ei,g :~ { (a cosh /i cos I/, a sinh /i sin z/) : < fi < fiQ, |i/| < i^o or Jtt — z/| < fo} 

D {(xi, 3:2) G -E : |a;i| > ai cosi^o}- 

Let 



1/2 



Pi^o ("^> '^) 



j j^y'^m.n) ^•^ 



Our particular interest in this appendix is in families of eigenfunctions that are ex- 
ponentially localised around the periodic orbit {(0,X2) : |a;2| < 02}. In particular we 
will show below that the family m„i^oj tj = 0, 1, ... is so localised; precisely, we will 
show that, for all vq g (0,7r/2), there exists /3 > such that ^2^^(777,, 0) = 0(e~^'^'") as 
m — )• 00. 

Noting (jA.5[) . we see that 

/ (um,nf dx ^ Aa^ j j (sinhV + sin^i^) (^Mc^^^(/x, (3'm,„)ce„(i/, (?,„,«) j d^idv. 
Thus, defining 



(\ 2 r^ 

Mc^^\^jL,qra,n)] d^, Is{m,n):^j (ce„(z^, g™,„))^ dz/. 



lO ^ ' Jo 

it holds that 

2/ (Ml -|-Mosin^i/o)-^i/o("^i"-) 
[p^g[m,n)) < 7-;-^^ , ^ , .^ .,_ 2 



Mil^/2{m, n) + Mo sin i^o(-?^77/2(w, ") - -^i^o ("^> ")) 

~ Iw/2{m,n) -/^„(TO,n)' 
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It is sufficient for the needs of this paper to estimate the asymptotics as ?7i — >■ oo 
of pjjg {m, n) for n = 0, and so we wih restrict our attention to this case. For this 
purpose, and abbreviating qmfl as q™ and fcm.o as fc^, recaU that ceo(i^, Qm) satisfies 
(jA.Sp with q = Qm, and with a = ao{qm,)- Now the asymptotics of the eigenvalue ao{q) 
as q ^>- oo are known. From 14 we have that 

a„{q) = -2q + q'^^ + 0{1) (A.9) 

as q -> oo. Thus we see that, for m large and with N = ceo(-,gm), the coefficient 
ao{qm) — 2g„jCOs2z/ of N{i') in (JA.3P is negative except in small neighbourhoods 
of ±7r/2 of length 0{qm ) = 0{km )■ It is this which causes the exponential 
localisation of Um.o around the periodic orbit. 

To see this localization completely explicitly, we will use the following lemma 
which depends on standard weighted space arguments (cf. [5] , [THl §3] ) . In this lemma 
and subsequently BC(R) denotes the set of functions : R — >■ R that are bounded 
and continuous and SC'=(R) := {(^ S BC{R) : (f)^^') S BC{M.) for j = 0, 1, . . . , k}. As 
usual, iJ*(R), for s > 0, denotes the standard Sobolev space of order s (which in this 
appendix we take to be a space of real- valued functions). For /3 > and a > let 
W/3,Q : R — > R denote the weight function 

W/3,q(s) = exp (— /3v a^ + 5^1, 

and, for 0G BC{R), let 

N 1/2 

{■Wl3^a{s)4'{s)fds} < OO. 



Lemma A.l. Suppose that p e i?C(R) and that, for some c> 0, p{s) < — c^ < 0, 
s e R. Suppose also that g G BC{R), v e BC'^(R), and 

v"{s)+p{s)v{s)= g{s), seR. 

Then, for < /3 < c and a > 0, 

\\vh,o.<^^\\9h,o.. (A.IO) 

Proof Suppose ^ > and a > 0. Define V € BC2(R) n ij2(M) by ■0(s) = 
^0,a('S)w(s), s € R. Then it is an easy calculation, abbreviating w^a as w, that 

^" -2—i,' +\p + 2[ — ] \ik^wg. 

w \ \w J w I 

Multiplying by a test function (/) and integrating by parts, we see that 

a(7A,<^)=6(<^), (f>eH'{R), (A.ll) 

where the bilinear form a : iJ^(R) x H^{M.) — > R and the bounded linear functional 
& : 7Ji(R) -j> R are defined by 

a(0, V) :- y" U'cb' + 2^i;'4> - Ip + 2 (^\ _ ^ ) ^0 ) d^, 

/OO 
wg(j)ds. 
-OO 
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Since w' /w, w" /w G BC{M.), the bilinear form a is bounded. For (f) e if^(R), 



2 



a{^,^)^l (0')2 + ^(02y-U + 2(^j -r_)0^)d. 



A'^2 



w 



2N 



where the last step follows by integration by parts, on noting that {w' /w)' = w" /w - 
{w' /w)'^. Since —p— {w' /w)^ > c^ — P^, a is coercive if /3 < c, with 

/ poo \ 1/2 

a(0, 0) > Ml where ||0||i :- (^y ((0')' + (c' " /?') 0') ^^ 
Applying the Lax-Milgram lemma, it follows from (lA.lip that 



V^^^^lkll/3,.<||^||l<||6||< 



I.9II/3,, 



y/c^^J^ 



where ||&|| denotes the norm of the linear functional b : H^{M.) — > M, with iJ-'^(R) given 
the norm || • ||i. Hence (jA.lOp holds for /3 > and a > 0, and so also for a = by 
the dominated convergence theorem. D 

To apply this lemma to (|A.3|) . let v := ceo(-, q) and define p S BCiM) by p(i^) = 
ao((j) — 2(7 cos 2i/, z/ € M. For c > write p a.s p = p^ + pt where p~ := min(— c^,p) 
and p+ :^ p — p~ , and set gc := ^ P^i'- Then w" +Pfrw = 5c, and applying Lemma 
lA. II with a = it follows that, for i/q e (0,7r/2), 

for < /3 < c. For — 2g < ao{q) + c^ < 2q we see that gdi^) = if |z^ — jtt\ < i^c, for 
some j e Z, where z^c G (0, 7''/2) is given by 

Thus, and since < p+ < ao(g) + 2q + c^ < Aq, it follows that 

\Ml,^2 re''^^igds)fds 
Jo 



.■_n J -i-lT + Vr 






g-2/3iy<: /.IT 

< 8(7^ ^:j^ / v^ ds. 





So 



^a ~| 1/2 g-/^(,.,-,.o) f f-^f^ 

v^ds\ <2y/2'q TbW v^ds 

J (c2 - /32) (1 - e-2/3^)l/2 [7o 



1/2 
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and choosing [3 — (?{vc — va)/{^ + ■\/l + {vc — v^Yc^ ) e (0,c), which minimises 

where 7 := (5/(l + \/l + <^), with 5 := c^{vc — vaY, provided that — 2q < ao(g) + c^ < 2g 
and < t'o < i^c- Thus, assuming that —ao{q)/{2q) e (—1,0) (which is certainly the 
case for aU sufficiently large q by (JA.9P ). and provided 

0<.o<icos-(^), (A.14) 

(|A.13P holds for < c < y/2q cos 2f — oq (q) . It particular, choosing 



c := -V29 cos 21/0 - aoiq) , (A.15) 

we obtain the following result in which the formula for 6 follows from S = c?[vc — 
i/q)^, with c given by (|A.15|) and Vc by (|A.12|) . recaUing that cos~^(— 1 + q) = tt — 
2 sin"^ v/c^' for < a < 2. 

Theorem A. 2. Let v := ceo{-,q) and suppose that — ao(g)/(2g) e (—1,0) (which 
certainly holds for all sufficiently large q). Then, provided vq satisfies (jA.14p . it holds 
that 

J 7(l-e-47-)i/2 Wg j 



where 4>q — it /2 — I'o, 7 = 5 /{I + VT+J), 

(5 = 9(sin2 00 - r(g)) (^^_^:,^-i f -y^sin^ 0o + 3r(g) j j . 

andr((7) := (1 + ao((z)/(2(7))/2. 

By (|A.9[) . r'(q) ~ g^^/^/2 as g -> 00. Thus, in the above theorem, as g — >■ cx) 
5 has the asymptotic behaviour 5 ^ qsin 0o (0o — sin~ (isin^o)) . Thus, and 
since sin~^ (isin^o) < \'t>Oi it holds that 5 > jq(j)Qsm'^ (po for all sufficiently large 
q. Further, 7 = (5^/^ — 1 + 0{S^^^^) as (5 — > 00. Thus the above theorem has the 
following corollary. 

Corollary A. 3. Let v := ceo(-,g). Then, for every vq € (0,7r/2), it holds for 
all sufficiently large q that 

"" 2,1'^' 2^2600 / 1 ^. ■ ,\\ r'^ 2,V^ 

V ds> < . — exp I --^0osm(/)o j < v ds "> , 

where (J)q — tt/2 — vq. 

Applying this corollary with q = qm = \{kma)'^, and recalling the bound (jA.8[) . 
and that a = aie, where e is the eccentricity of the ellipse, and that (Jjq/ sin c/jq < 7r/2, 
we see that, for every I'o G (0, 7r/2), it holds for all sufficiently large m that 

Pyo(m, 0) < V2 7re exp f --efcmai(/)o sin^o ) , (A.16) 

where c/jq = "k 12 — vq . 
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